#! /bin/bash
#
# Generates a set of circuit-defined rotations.
#
# Get started
#echo "  "
#echo "-------------morafr circuit closer--------------"
#echo "  "


echo ">" > grneur_temp.dat
echo ">" > eurnam_temp.dat

d2u grneur_temp.dat
d2u eurnam_temp.dat

GRNEURINVLIM=`cat C:/GMT/2019/GRNEUR3/inversion/data/iter.dat | awk '{printf "%02d",$4}'`
EURNAMINVLIM=`cat C:/GMT/2019/EURNAM/inversion/data/iter.dat | awk '{printf "%02d",$4}'`


echo "Using these rotations:"
echo C:/GMT/2019/GRNEUR3/inversion/data/poles_$GRNEURINVLIM.dat
cat C:/GMT/2019/GRNEUR3/inversion/data/poles_$GRNEURINVLIM.dat
echo "..... and ....."
echo C:/GMT/2019/EURNAM/inversion/data/poles_$EURNAMINVLIM.dat
cat C:/GMT/2019/EURNAM/inversion/data/poles_$EURNAMINVLIM.dat

echo " "

LISTLENGTH=`awk '{if (NR==1) print $2}' C:/GMT/2019/GRNEUR3/inversion/data/poles_$GRNEURINVLIM.dat`
COUNTLIM=`echo $LISTLENGTH "1" | awk '{print $1-$2}' | head`
#
#echo $GRNEURINVLIM $EURNAMINVLIM $LISTLENGTH $COUNTLIM
# generate the list of rotations for which the addition will be done
#
i=0
while [ $i -le $LISTLENGTH ]
do
 CHRONLIST[$i]=`awk '{if (NR==i+2) print $4}' i="$i" C:/GMT/2019/GRNEUR3/inversion/data/poles_$GRNEURINVLIM.dat`
 i=$(( $i + 1 ))
done
#
echo "Will complete the circuit for these chrons: "${CHRONLIST[*]}
echo ">" > chroncol.dat

d2u chroncol.dat


# while loop
i=0
while [ $i -le $COUNTLIM ]
do
  CHRON=`echo ${CHRONLIST[$i]}`

  awk '{if ($4==CHRON) printf "%10.3f %10.3f %10.3f\n", $1,$2,$3}' CHRON="$CHRON" C:/GMT/2019/GRNEUR3/inversion/data/poles_$GRNEURINVLIM.dat >> grneur_temp.dat
  awk '{if ($4==CHRON) printf "%10.3f %10.3f %10.3f\n",$1,$2,$3}' CHRON="$CHRON" C:/GMT/2019/EURNAM/inversion/data/poles_$EURNAMINVLIM.dat >> eurnam_temp.dat
  echo $CHRON >> chroncol.dat
  
 i=$(( $i + 1 ))

done

tail -$LISTLENGTH grneur_temp.dat > crap_a
tail -$LISTLENGTH eurnam_temp.dat > crap_b
#tail -$LISTLENGTH indant_temp.dat > crap_c

d2u crap_a crap_b

paste crap_a crap_b > circstep1_add.dat
d2u circstep1_add.dat
#rm crap_a crap_b

#adds the rotations together

echo "Expecting you to have sent me this format for adding rotations together: "
echo " "
echo "LATITUDE LONGITUDE ANGLE LATITUDE LONGITUDE ANGLE"
echo " "

# Reformat for the fortran
#
awk '{printf "%13.6f%13.6f%13.6f%13.6f%13.6f%13.6f\n",$1,$2,$3,$4,$5,$6}' circstep1_add.dat > temppol.dat

# Now the fortran
####################################
./poles_add temppol.dat tempout.dat
tail -$LISTLENGTH tempout.dat > crap_c

#paste crap_d crap_c > circstep2_add.dat
#awk '{printf "%13.6f%13.6f%13.6f%13.6f%13.6f%13.6f\n",$2,$1,$3,$4,$5,$6}' circstep2_add.dat > temppol.dat
#poles_add temppol.dat tempout.dat

echo "Written these rotations:"

paste tempout.dat chroncol.dat | awk '{if (NR==1) print "grnnam "LISTLENGTH; else print $2,$1,$3,$4}' LISTLENGTH="$LISTLENGTH" > ../grnnam_pol_circ.dat

#echo "... to /Users/geagles/ausant_circuit/ausant_pol_circ.dat"
echo "..."
echo "Bonzer"

